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ABSTRACT 

We consider the deformation that has recently been observed in the inner part of the 
circumstellar disk around [3 Pictoris with the HST. Our recent ground based adaptive 
optics coronographic observations confirm that the inner disk is warped. We investi- 
gate the hypothesis that a yet undetected planet is responsible for the observed warp, 
through simulations of the effect of the gravitational perturbation due to a massive 
companion on the disk. The physical processes assumed in the simulations are dis- 
cussed: since the observed particles do not survive collisions, the apparent disk shape 
is driven by the underlying collisionless parent population. The resulting possible pa- 
rameters for the planet that are consistent with the observed disk deformation are 
reviewed. 

Key words: stars: individual: /3 Pictoris - circumstellar matter - planetary systems. 



1 INTRODUCTION 

An important discovery by IRAS was the detection of IR 
excess due to cold material around a number of Main Se- 
quence Stars. Among them, f3 Pictoris exhibits one of the 
largest excess (Aumann, 1985). Shortly after the IRAS re- 
sults, a disk of dust was imaged around this star (Smith 
and Terrile, 1984). Since then, the disk has been extensively 
studied , as it is expected to be related to planetary sys- 
tems, possibly in a state of evolution different from that of 
our Solar System. 

The disk composition is complex. Micron sized grains 
are detected through scattered light (Kalas and Jewitt, 1995 
and ref. therein; Mouillet et al, 1996) and thermal emission 
(Lagage and Panin, 1994). The presence of small amounts of 
submicron grains is inferred from 10 /im spectrophotometry 
(Knacke et al, 1993). Larger grains (mm sized) are detected 
through photometry at mm wavelengths (Chini et al, 1991; 
Zuckerman and Becklin; 1993). Kilometer-sized bodies have 
also been proposed to account for the very peculiar spectro- 
scopic variability of /3 Pictoris (Lagrange et al, 1987). Gas is 
detected spectroscopically through the presence of absorp- 
tion lines. However the gas to dust ratio is probably < 1., 
ie much less than that in the environment of young stars 



such as T Tauri, or in star forming regions. Most of the gas 
is probably confined very close to the disk (Lagrange, 1995 
and ref. therein). 

The total mass of the disk is very uncertain as most of 
it comes from the largest bodies. About one Earth Mass is 
necessary to account for optical to mm observations. If the 
particle size distribution follows a law of the form, 



dn(a) 



3 da, 



with n(a) being the number of particles of size greater than 
a, up to kilometer size bodies, the mass would be several 
tens of Earth Masses (Backman and Paresce, 1993). 

The origin and evolution of this system is still the sub- 
ject of active research and one point of investigation is the 
possible presence of planets within the disk. From the theo- 
retical point of view, planets are expected to be formed out 
of the the circumstellar disk that accompanies the star dur- 
ing the process of formation. The timescale estimated for 
this process (see Lin and Papaloizou, 1985) is characteristi- 
cally less than the estimated lifetime of /3 Pictoris ~ 2.10 8 y 
(Paresce, 1991). Thus is is reasonable to suppose that plan- 
ets have had time enough to form. 

Obtaining direct evidence for the existence of plan- 
ets around early type stars such as /3 Pictoris is far from 
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straightforward. Direct imaging is still beyond current obser- 
vational capability because of the high contrast between the 
star and the planet. Photometric variations have been ob- 
served (Lecavelier et al, 1995) which are consistent with the 
occultation of /3 Pictoris by an orbiting Jupiter-like planet. 
But this has not been firmly established (observation of an- 
other occultation is needed to sweep all doubts away). Ra- 
dial velocity studies, such as those performed to detect giant 
planets around solar type stars such as 51 Peg, 47 UMa, and 
others (Mayor and Queloz, 1995, Marcy and Butler, 1996) 
would not provide enough accuracy to detect giant planets 
at distances of about a few AU from these ususally rapidly 
rotating early type stars. 

Signatures of planets can result from their gravitational 
effect on the disk. The observed clearing of the inner region 
of the disk has been attributed to the gravitational effect of 
an orbiting planet (Roques et al, 1994). But other physical 
processes can explain the observed dust distribution as well. 

The gravitational influence of a planet has been invoked 
in order to explain the observed high rate of cometary in- 
fall (about 1000 per year corresponding to ~ 10~ lr> Moy^ 1 ). 
Beust and Morbidelli (1995) propose that mean motion res- 
onances between the orbits of kilometer sized bodies and a 
planet in a longer period eccentric orbit can produce the 
right quantity, and orbits of infalling bodies. They predict 
long term variabilities which are currently being tested via 
a spectroscopic survey. Such a model requires a steady flux 
of kilometer sized bodies from the outer to inner regions of 
the disk to occur as a consequence of collisions. 

Finally, HST data provide high S/N data on the light 
scattered from the the inner part of the disk down to 25 
AU from the star (Burrows et al, 1995). They revealed a 
slight inclination (~ 3 degrees) of the inner disk midplane 
up to about 50 AU, to the outer disk midplane. Burrows et 
al (1995) proposed that a planet on an orbit inclined to the 
outer disk plane might be the most probable explanation for 
this observation. In section 2, we show recent ground-based 
observations imaging the same region of the disk for com- 
parison. In section 3 we review the physical properties of the 
disk to determine how an hypothetical planet would act on 
it. This leads to 3-D numerical simulations of the simultane- 
ous behaviour of large numbers of test particles represent- 
ing a disk population of kilometer sized objects under the 
gravitational field of the star and massive companion on an 
inclined, and possibly eccentric orbit. In section 4 we present 
the results. We finally discuss the possible parameters for the 
companion consistent with its presence being able to explain 
the brightness asymmetry of the inner /3 Pictoris disk. 



2 OBSERVED ASYMMETRY OF THE INNER 
PART OF THE j3 Pictoris DISK 

2.1 Adaptive optics observations of the inner disk 

No warp has been detected in the classical coronographic ob- 
servations outside 80 AU. In order to observe the disk closer 
to the star, one needs high angular resolution observations, 
ie either space observations or atmospheric effect corrected 
images. Recently, adaptive optics coronographic images of 
the (3 Pictoris disk enabled us to detect the disk in the near 
infrared (NIR) through scattered light down to about 25 
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Figure 2. The observed vertical deformation of the disk above 
the outer disk midplane. The axes are marked in AU. The verti- 
cal deformation is measured from the centroid of the brightness 
distribution in the direction normal to the outer disk midplane. 
The value zero is adopted when the measurement is not possible 
(closer than 25 AU and further than 110 AU). The solid line is the 
linear fit for the inner part. We define the extent of the deforma- 
tion as the distance at which the slope of the fit starts to decrease. 
The extent of the deformation is 50 AU, and the corresponding 
inclination 3 degrees 



AU from the star (Mouillet et al, 1996). Similar resolution 
was obtained with HST in the optical range (Burrows et al, 
1995). 

We performed new observations of the f3 Pictoris disk 
with a coronograph coupled to the ESO adaptive optics sys- 
tem ADONIS on January 5, 1996 in the J band (Fig. 1). The 
inner region (30-80 AU) is observed with an angular reso- 
lution similar to that of HST. However, the correction for 
atmospheric turbulence makes it more difficult to subtract 
perfectly the stellar light which remains after the mask, so 
that some residuals of the PSF temporal variations remain 
on the image. Details of the observing method and of the 
reduction procedure can be found in Beuzit et al (1996). 

From Fig. 2 we see that the intensity peak lies on one of 
two lines according to whether it is nearer or further than 50 
AU from the centre. Thus the disk midplane inside 50 AU 
is inclined at about 3 degrees to the outer disk midplane. 
The effect of the warp is hardly seen, but nonetheless it is 
consistent with the HST data in terms of amplitude and 
extent. 

These observations confirm previous results (Mouillet 
et al, 1996; Burrows et al, 1996): the surface brightness pro- 
file I(r), which is measured to be very steep far from the 
star (Kalas and Jewitt, 1995), gets flatter and flatter as the 
distance to the star, r, decreases from 100 AU down to 25 
AU. For r > 100 AU, I(r)oc r~ 3 ' 5 . For 50 AU < r < 100 AU, 
I(r)oc r -1 ' 3 . For r < 50 AU, the profile is even flatter. (Fig. 
3). Consequently, the form of the optical optical depth as a 
function of r, given by r(r) = 5 10~ 3 (r/100) -1 ' 7 , which is 
valid for r > 100 AU, cannot be used for smaller values of r. 
The inferred form depends on the precise model involving a 
combination of thermal IR and scattering observations, but 
it converges on a progressive flattening up to a maximum 
value of 10~ 3 -10~ 2 at around r — 40 AU, thereafter de- 
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Figure 1. Adaptive optics observations in the J band of the inner disk of fi Pictoris, January 05 1996. The axes are marked in AU. 
Inward of 50 AU, the disk midplane is inclined with respect to the outer disk midplane (dashed line) 
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Figure 3. Midplane surface brightness of the disk 

creasing inwards. Very little residual dust is expected closer 
than 15 AU, since it would induce IR excesses in NIR which 
are not observed. 

2.2 The planetary hypothesis 

The deformation of the inner disk midplane relative to the 
outer disk midplane is not the only asymmetry observed 
in the j3 Pictoris disk (Kalas and Jewitt 1995; Lagage and 
Pantin, 1994) . All these asymmetries are difficult to explain 
in the context of a low mass disk consisting of many small 
objects under a dominant potential arising from a central 
point mass, since any azimuthal inhomogeneities would be 
expected to disappear on an orbital timescale (a few 1000 
years at 100 AU). A single physical process is unlikely to 
be responsible for all of them since very different spatial 
scales are involved together with different geometries. For 
instance, anisotropic scattering properties together with a 
slightly inclined disk were proposed to explain the so called 
"butterfly" asymmetry (Kalas and Jewitt, 1995: the mid- 
plane is not an axis of symmetry of the brightness distribu- 
tion). But such an explanation is unable to count for any 
radial-dependent asymmetry. Another idea is that the local 



production of very small grains as collision products followed 
by their efficient accelaration by radiation pressure could 
produce catastrophic blow out in some random part of the 
disk (Artymowicz, 1996). Such a mechanism could explain 
any kind of asymmetry without the need of a planet. But the 
strong antisymmetry ( global azimuthal structure) of the ob- 
served warping with respect to the central star make this an 
implausible explanation for it. In this paper we investigate 
the possibility that the warping asymmetry is produced by 
an as yet undetected planet in order to see what constraints 
might be put on its parameters. Although we are not able to 
claim this as a unique explanation, we are able to show that 
the required parameters are reasonable and not at present 
in conflict with constraints arising from the lack of observa- 
tional detection through either imaging or spectroscopy. 



3 NUMERICAL SIMULATIONS OF THE 
INFLUENCE OF A PLANET 

3.1 Physical processes in the (3 Pictoris disk 

The dynamics in the disk is dominated by the gravitational 
force of the central star to which may be added effects due to 
any possible perturbing planets. The mass of observed disk 
material is negligible. Radiation pressure is important for 
small grains less than a characteristic size a c = (Arty- 
mowicz, 1988). For such grains the force due to radiation 
pressure exceeds gravity so that they are on unbound hy- 
perbolic orbits. Accordingly their lifetime is on the order 
of the crossing time of the disk and thus very short. The 
Poynting Robertson effect is negligible as is gas drag due to 
the small quantity of gas present (Artymowicz, 1995). 

The main quantity determining the grain collision fre- 
quency is the optical depth of the dust, r. The mean time 
between collisions is given by 



Poi 



tc = 



where Porb is the local orbital period. For the 
/3 Pictoris disk, the optical depth has been estimated to 
be typically 5 x 10~ 3 for small size scattering particles. 
Supposing that the particle size distribution is given by 
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dn(a) oc a~ 3,5 da, ( which might arise from collision pro- 
cesses) the smallest micron sized grains contribute most of 
the optical depth. The collision time at WOAU is then typ- 
ically W 5 y for such grains. The characteristic impact veloc- 
ities are ~ v r = iv or b ~ lfcm/s, where v or b is the orbital 
velocity and i is a characteristic inclination in radians as- 
sociated with a grain orbit. The expected impact velocities 
are sufficient to destroy the colliding particles. Therefore the 
lifetime of the smallest particles in bound orbits is on the or- 
der of the collision time ~ 10 y at 100 AU and even shorter 
closer than 100 AU, in any case much less than the age of 
the system. Thus they need to be replenished by the destruc- 
tion of larger particles. For normal solid state densities, the 
surface density of disk material contained in small particles 
is E ~ a c r, which implies a mass of ~ 10~ 9 M© contained 
within WOAU, and a mass loss rate of around 10 _14 Mqj/ _1 . 

Assuming that collisions will determine both the par- 
ticle size distribution through fragmentation processes 
(dn(a) oc da) as well as the particle lifetimes, most 

of the mass in the system will reside in the largest bodies 
which also contribute least to the optical depth. It is natural 
to regard the larger bodies, as a primary source of material 
to drive the global mass loss rate of the system. After a (de- 
structive) collision, the particle ejection velocities, as viewed 
in the centre of mass frame, are negligible compared to the 
orbital velocity of the parent particles. Consequently, the 
global distribution of any given size particles is not driven 
by interaction processes but is closely related to that of the 
much bigger parent bodies. This means that the disk shape 
of the small particles detected in scattered light is driven by 
the kinematic effect of gravitation on long-lived (old as the 
system) particles . 

Because of the detection of this disk, the collision time 
of the biggest bodies is expected to be larger than the age 
of the system. As an illustrative example, if the largest 
sized bodies are km sized, the given size distribution in- 
dicates that ~ 3.1O _5 M0 should be contained in these bod- 
ies and the characteristic collision time between them to be 
~ 3.10 9 y. If a significant amount of the material involved 
in collisions ends up in small particles, then the supply rate 
could be as much as ~ 10 _14 M©j/ _1 . Note too that collisions 
between the km sized bodies also result in an inward radial 
migration of the distribution much as for viscous accretion 
disks or planetary rings. The effective diffusion coefficient 
is D c ~ H 2 /t c , H being the vertical thickness and t c being 
the collision time. For H = O.lr, and t c = W 9 y, this gives 
a diffusion timescale of r 2 /D c ~ W n y, corresponding to 
a mass flow rate of ~ 3.10 - 16 Mgy -1 . Apart from lifetime 
considerations, a large number of such km sized bodies is 
expected in the /3 Pictoris disk to explain the spectroscopic 
detection of highly redshifted variable gas close to the star 
(Beust et al. 1996). Were even larger bodies to exist around 
the star, like the kilometer sized objects they might also pro- 
vide a source of small particles, but they would be expected 
to have a similar or perhaps even longer effective collision 
frequency, and thus behave in a similar way from the point 
of view of our analysis. 

Accordingly we model the disk as a system containing 
a large number of long-lived objects, which as far as the ma- 
jority is concerned, are collisionless. We then suppose that, 
as they do not survive more than one collision,the small 
particles reflect the distribution of the large bodies. We are 



Table 1. Parameter sets of the simulations in model units (see 
text) 



Model 


M/M* 


D 
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Run time 


01 


lc-3 


1. 


0. 
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le5 


02 


4e-3 


0.5 


0. 


3 


le5 


03 


9c-3 


0.33 


0. 


3 


3e5 


riA 
u^t 


OC-O 


1 

1 . 


n 
U. 


Q 
O 


ICO 


05 


lc-2 


0.7 


0. 


3 


le5 


06 


lc-2 


1. 


0. 


3 


le5 


07 


lc-3 


0.7 


0.1 


3 


7.5c5 


08 


lc-3 


0.7 


0.3 


3 


3e5 


09 


lc-3 


0.7 


0.5 


3 


3e5 


10 


lc-2 


0.7 


0.1 


3 


3e5 


11 


lc-2 


0.7 


0.3 


3 


3e5 


12 


lc-2 


0.7 


0.5 


3 


3e5 


12b 


lc-2 


0.7 


0.5 


6 


3e5 



thus concerned with kinematic patterns produced by non- 
interacting particles. 

3.2 Numerical simulations 

The numerical simulations are based on a collisionless adap- 
tation of the SPH code described in Larwood and Papaloizou 
(1997, and references therein). In this version pressure and 
viscous forces are removed so that we consider a purely 
kinematic model. The code follows the velocity and posi- 
tion coordinates of typically 15,000 particles over thousands 
of revolutions of the perturbing companion. The particles 
were initially set up in circular motion in a disk configura- 
tion with aspect ratio H/r = 0.1. They were inserted, using 
a random number generator according to the following E 
distribution (in model units) to schematically represent the 
matter distribution in the (3 Pictoris disk (since the particles 
are collisionless, the precise form of E(r) has no effect on the 
amplitude and extension of the consequent warp) : 

- for 1.2 < r < 4.2, E(r) ex r~°' 5 

- for 4.2 < r < 6, E(r) cc r" 1B . 

The model unit of length corresponding to r = 1, R u „u, is 
arbitrary, whereas the corresponding time unit t un it is the 
inverse orbital frequency at Runit- 

The companion is initiated in a generally eccentric or- 
bit such that, in a Cartesian coordinate system with origin 
at the center of mass, and (x,y) plane coinciding with the 
initial disk midplane, both the planet and the line of nodes 
are on the x axis at time t = 0. Parameters associated with 
the simulations are, the mass ratio of the planet and central 
star, M/Mt, the semimajor axis, D, eccentricity, e, and incli- 
nation of the planet's orbit, i, to the disk midplane. Finally 
the run time since initiation with the orbital configuration 
as described above is also an important parameter as this 
determines the extent of the warped deformation. Table 1 
gives the chosen parameters for various simulations. 

In order to compare the simulation results with obser- 
vation, we have derived the apparent surface brightness dis- 
tribution associated with the simulation particles. That is 
the integration along a line of sight of the scattered stellar 
flux from each particle (Fig. 4). We note that a companion 
in an inclined orbit breaks the system's azimuthal symme- 
try. The orientation of the line of sight to line of nodes is 
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Figure 4. Apparent surface brightness distribution associated with the simulation particles. The represented data set is model 05 at 
run time = 10 5 . For each projection plane, the axes are scaled in model units 



unknown in the case of f3 Pictoris, and should be considered 
as a free parameter in fitting the observations. 



4 RESULTS 

4.1 Simple description of the effect of a 

companion on the vertical elevation of the 
disk 

Adopting the cylindrical coordinates (r,<p,z), equivalent to 
our Cartesian system, the vertical elevation of a razor thin 
disk moving under the influence of a companion in an in- 
clined circular orbit with small i satisfies 



orbits, tp = Qt + const., and for small e and i 



d 2 , r>2 



3 GMD 2 



sin(2i) sin(y). 



(1) 



Here fi z is the frequency of vertical oscillations and the term 
on the right hand side is the secular vertical acceleration due 
to the companion. This has been expanded to lowest order 
in D/r, which becomes an increasingly accurate procedure 
further out in the disk (Larwood and Papaloizou, 1997). The 



precession frequency lo p — fi — fi z , is such that 



< fi. 



where fi is the angular velocity. For particles in near circular 



3 GMD 
w K = — — - — . 

p 4 fir 5 

Using the above in (Q), we obtain the approximate so- 
lution with zero elevation at t — 0, 

3 GMD 2 sin(2z) 
2 = ~b o 3 (sm(vj) - sin(¥3 - Wpt)) . (2) 

In this purely kinematic model, the inclined orbit of 
the companion breaks the midplane symmetry and forces 
the test particle orbits to precess around the companion's 
orbital angular momentum axis. If phase mixing due to the 
r dependence of lo v is complete (\u} p t\ 3> 1), this process 
thickens the disk so that the aspect ratio appears to be 
equal to the orbital inclination. The vertical amplitude of 
the deformation is thus directly related to the inclination of 
the companion's orbit. 

This relation is verified numerically for the simulation 
data (Fig. 5). In addition, we note that the apparent ampli- 
tude of the warp depends on the orientation of the line of 
sight (Fig. 6). The observability of the warp is maximized 
when it (assumed to lie in the outer disk midplane) is at 
about 30 degrees to the line of nodes. 
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Figure 5. Linear dependence of the amplitude of the warp on 
the inclination of the companion's orbit. The axes are scaled in 
model units. The amplitude of the warp is quantified in the same 
way as for observational data (see Fig. 2). The effect is twice as 
large for an inclination of 6 degrees (dashed line, from model 12b), 
compared to an inclination of 3 degrees (solid line, from model 
12), all other parameters equal 
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Figure 7. Radial extent of the deformation. The extent of the 
warp (in model units) is measured in the same way as for the 
observational data (Fig. 2). It is displayed versus MD 2 t in a log- 
log plot. The simulation data comes from zero eccentricity models 
(modelOl: +, model02: *, model04: O, model05: A, model06: □) 
and a line of sight orientation of 60 degrees (where the warp 
observability is unambiguous). Similar results arc obtained with 
other models and orientations. The best fit power law index is 
numerically measured to be 0.29, in good agreement with the 
theoretical expectation (2/7) 



that we expect the scaling relation 

r cc (MD 2 t) 2/r . 

This form occurs because the propagation rate is driven by 
the tidal force which is proportional to AID 2 . 

The radial dependence of the precession frequency, 
namely 



lu v cc MD r 



-7/2 



has been checked in the simulations (Fig. 7) 



-3 F 1 

-100 100 

line of sight orientation (deg) 

Figure 6. Observability of a warp as a function of the line of 
sight orientation. The x-axis gives the orientation of the assumed 
line of sight to line of nodes. The y-axis gives the corresponding 
apparent inclination of the warp in simulated data (with a 3 de- 
grees inclined companion orbit), measured in the same way as for 
the observational data (Fig. 2). 



4.2 Radial extent of the deformation 

A deformation corresponding to a thickenned disk is estab- 
lished very fast close to the star and then propagates out- 
wards. As indicated above by (^) the location of the end of 
the deformation very roughly corresponds to the place where 
the inverse precession frequency is roughly equal to the run 
time (\uJ p t\ ~ f). But note that this will overestimate the 
location of the edge of the deformation because many pre- 
cession periods are needed to ensure enough phase mixing 
for the full deformation to be produced. But the number of 
precession periods required should not depend on radius so 



4.3 Companion eccentricity and particle inflow 
towards the star 

We now explore the consequence of a possible eccentricity 
of the companion orbit. Increasing eccentricities induce very 
little direct observable effect for a disk seen edge-on (Fig. 
8,9). Consequently, this value is poorly constrained by the 
observations of the /3Pictoris disk. Yet, the azimuthal distri- 
bution of particles is affected, with the formation of spiral- 
like features. The stronger the eccentricity, the sharper such 
features get. 

We also notice that the eccentricity of some disk parti- 
cles gets increased to high values with the result that they 
are thrown towards the star and end up inside the compan- 
ion's orbit. Out of the initial number of 15354 particles with 
r > 1.2, a few get inside the companion orbit (D — 0.7) in 
model 11 (e = 0.3), and a few dozen in model 12 (e = 0.5). 



5 POSSIBLE PARAMETERS FOR THE 
PERTURBER 

The comparison between the simulation results and observa- 
tions allows derivation of possible parameters for the com- 
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Figure 8. Same as Fig. 4, for modcllO at run time 10 J : the companion orbit eccentricity is 0.1 



panion orbit. First, the observed inclination of 3 degrees 
requires an orbit inclination of 3-5 degrees to the outer 
disk midplane, depending on the angle between the line of 
sight and the line of nodes ( these assumed to be coplanar) . 
This angle is poorly constrained but should not lie in the 
[—80°,— 50°] (modulo 180°) range so that the deformation 
is observable. 

Second, the companion orbit is not well constrained by 
the observation of an edge-on disk, so that it can be consid- 
ered as a free parameter in the simulations. This could have 
an important role in feeding material into the inner disk. 
This, taken together with the general inward migration of 
the km sized bodies may be related to the cometary infalls 
observed via spectroscopy, at the rate of about 1000 km- 
sized bodies per year. Gravitational perturbations are likely 
to be involved in these processes. A massive companion in- 
voked to explain the imaging data, could then also provide 
at least part of an explanation for these infalls, as long as 
its eccentricity is significantly larger than 0.1. We also re- 
mark that a moderately large eccentricity is required if an 
eccentricity pumping mechanism, through mean motion res- 
onances, is supposed to work on particle orbits interior to 
the planet. 

Finally, observations show a radial extension of the warp 



of 50 AU. If we arbitrarily choose the model unit R U nit = 
10 AU, then the corresponding model time unit is : 



G 



1/2 



5-2 y. 



According to the simulations (Fig. 7), the observed extent of 
the warp (log(R/ Runit) = 0.7) requires then the condition : 



M 



V Runit J 



t 



3000. 



tunit / tunit 

Assuming that the companion formation is rapid, so that 
the propagation time of the warp is close to the age of 
the system, t ~ 2 10 s y. Then a planet with MD 2 = 
10~ 4 M„(10AU) 2 can account for the observed warp (Fig. 
10). In the framework of these simulations, the compan- 
ion orbit is interior to the observed disk so that D < 20 
AU. The absence of radial velocity variations larger than 
about 1 km/s in our spectroscopic data gathered since 1984 
implies that M 2 /Dav < 5. 10~ 4 . Meanwhile, the non di- 
rect detection of any close by companion does not con- 
strain further the range of possible parameters. Finally, 
the possible parameters define a giant planet with a mass 
10~ 5 < M/M+ < 10~ 2 , located at a corresponding respec- 
tive distance from the star 204(7 < D < 1417. Such a 
planet may be detected in photometry if it crosses the line 
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Figure 9. Same as Fig. 4, for modcll2 at run time 
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Figure 10. Possible parameters M and D for the gravitational 
perturber in a log-log diagram. The constraint derived from the 
extension of the warp is represented by the dashed line. Shaded 
regions arc forbidden because of observational constraints (see 
text) and the assumption that the companion orbit is inside the 
disk 



of sight, which requires that the earth is in the plane of the 
planet orbit to within a precision of 0.06° /Dau- Actually, 
Lecavelier et al (1995) propose a Jupiter-like planet to ex- 
plain short term light variations, which is within the present 
possible parameter range. 

These quantitative results on the possible companion 
parameters directly depend on the system lifetime, t, which 
is still unprecisely estimated. However, Crifo et al (1997) 
derive from the stellar photometry and the new Hippar- 
cos data that the system may not be younger than W 7 y. 
In this case, our corresponding constraint would shift to 
MD 2 = 210~ 3 M»(10AU) 2 . This would restrict the pos- 
sible parameters to more massive companions (10 -3 ' 5 < 
M/M* < 10~ 2 ) between 20 AU and 3 AU. 



6 CONCLUSION 

From visible as well as NIR imaging observations, the disk 
around /3 Pictoris is detected through scattered light from 
/jm sized grains. These grains are short-lived because of de- 
structive collisions and the action of radiation pressure. As- 
suming that their distribution is the same as an underlying 
parent population of km-sized bodies, the apparent matter 
distribution within the disk can be derived from numerical 
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simulations of collisionless particles. Such simulations are 
able to reproduce the observed warp in the inner part of the 
disk as a result of the effect of gravitational perturbations 
due to a planet in an orbit inclined at 3-5 degrees to the 
outer disk midplane. Such a planet should be located be- 
tween 1 and 20 AU, with a corresponding mass respectively 
between 10~ 2 and 10~ 5 M*. 
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